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Abstract 

We study the cross section for dijet production via direct photons, as seen at the HERA 
collider. Rather than assuming the standard factorisation formula for the cross section 
(where the incoming partons are assumed to be on-shell), we use the /^-factorisation 
formula. This formalism ought to be more appropriate for jet production in regions where 
the cross section is sensitive to the small- a; parton distribution functions. We study the 
rapidity distribution and the azimuthal angular distribution of the jets and compare with 
the traditional approach. Our calculations reveal a breakdown of the traditional approach 
as one moves into the regime of small x. 



In a previous paper Q, we illustrated the importance of studying the cross section for dijet 
production via direct photons at the HERA collider in order to gain information regarding the 
gluon content of the proton at small- a;. Since that paper was published, the ZEUS collaboration 
has shown that it is possible to separate their events into 'direct photon' events and 'resolved 
photon' events using a cut on the fraction of the photon energy which enters the hard scatter 
0| . This encourages us to study further this important process and so, in this paper, we 
wish to return to the theoretical evaluation of the dijet cross section in the context of the 
'/^-factorisation' prescription, which ought to be used for small x ~ v\l s (Pr is a typical jet 
transverse momentum and s is the ep centre-of-mass energy). 

At large x ~ Pt/ s j we know that the lnpr terms, which are present in the perturbative 
expansion of the cross section, can become significant and need to be summed - leading to the 
violation of the scaling predictions of the parton model. The summation of the leading loga- 
rithms (which arise as a result of multi-parton branching before the hard scatter) is performed 
using the DGLAP (Dokshitzer, Gribov, Lipatov, Altarelli, Parisi) equations ||. The infra-red 
singularities (reflecting the presence of long distance physics) which arise in the calculation of 
these QCD corrections can be separated from the short distance contributions and thus we are 
led to cross sections which factorise (schematically) in the following way [Q: 



The function F contains the (unknown) long-distance physics whilst only short distance physics 
is contained in the 'hard' subprocess cross section a. It is usual to choose the factorisation 
scale (/i 2 ), to equal the hard scale (Q 2 ), in order to subsume the leading logarithms within the 
function F. Central to this factorisation formula is the assumption that the 'hard' cross section 
can be computed assuming that the incoming partons have negligible transverse momenta (and 
virtualities) compared with the scale of the hard process. We shall call this the 'Q 2 -factorisation' 
approach. 

For small enough x we expect that the In a; terms, which are present in the perturbative 
expansion, will become dominant. We refer to the summation of these leading logarithms in 
1/x as the BFKL (Balitsky, Fadin, Kuraev, Lipatov) summation ||. If the logarithms in x are 
indeed dominant, then we are no longer allowed to use the formula of eqn. ([]]). We must use 
the /cT-factorisation formula of refs.0, [7], ||, which states (again schematically) that 



F(...,/i 2 )®a(...,/i 2 ,g 2 ) 



(i) 



a ~ F(..., k 2 Ti , /i 2 ) <g> &(..., k 2 Ti , fj 2 , Q 2 ) 
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The scales, kTi, are the transverse momenta carried by the partons entering the hard scatter. 

In the specific case of direct photoproduction of dijets the (^-factorisation formula can be 
written: 

, d ? , 2 = z fi/e{z) xg(x,Q 2 ) da ^ 9 Z" q V (3) 
dy 1 dy 2 dpT dt 

The final state parton rapidities are denoted y\ and y 2 , and pt denotes their transverse momen- 
tum in the 7p centre-of-mass frame. The azimuthal angular distribution is trivial since the jets 
are produced back-to-back. The photon flux is determined by the function, f 7 / e (z). The long 
distance physics is contained in the gluon structure function, xg(x, Q 2 ), and the short distance 
physics is contained in the cross section: 

dajjg -> qq) = 2Tcaa s / £ u\ 
dt ~ s 2 \u t) 

For simplicity we show only the gluonic contribution to the cross section (it is this contribution 

which dominates for jets produced in the electron direction, i.e. small x). The subprocess 

Mandelstam variables are denoted s, t and u and satisfy the usual kinematic relations. 

The corresponding formula in the language of &/r-factorisation can be written thus ||: 

d 2 a dir . d 2 ^ i2^ d °{lQ*^W) 

—— — = / d P t <)(Pt + Pt - k Tj J(v,p T ,k T ,0) zf j/e (z) ®{x,k T ) - z 

dy\dy 2 J tx tx dt 

(5) 

The final state partons have transverse momentum vectors pt and p T which are no longer equal 
and opposite for non- vanishing gluon momentum, kx- We can use the delta function to perform 
the p T integrals. One of the two remaining angular integrals can also be performed (since the 
integrand only depends upon 6, the angle between kp and px). We can therefore write: 

pp = I dpidki f. j(»y T , 4, e) Z f 7/e (z) $(*, %) Ml9 *r m) (6) 

dy±dy 2 J 2tt dt 

For the 'hard scatter' cross section we use HI: 



da{^g* — > qq) 2naa s ft u 2k\s Ylk\s \ 

It = (I + k\Y + ~W + (u + £) 2 J () 

J is the Jacobean and v = exp (yi — y 2 ). Some important kinematic relations follow. 

Pt — Pt + &r — 2/c T p T cos 9 (8) 
The momentum fractions x and z can then be written 

x = ±-(p T e-y-+pi T e-y*) 

' (p T e yi +p' T e y2 ) (9) 



2E, 



and 

-i = p 2 T + PTP ' T /u (io) 

Also, s = AxzE e E p and s + i + u + k^ = 0. We can write the Jacobean explicitly as: 



n-l 2 .2 q\ _ PtPt 



' 2 A/ T A 2 \ A 2 

v 2 + v ^ + + — 

\Pt PtPt J Pt 



where A 2 = (p T + p? - k\) /2. 

To regain eq.([3J) from eq.(|S|), one assumes that the integrand can be approximated by its 
value at k\ = (the 9 integral is then trivial) and that the maximum value of k\ = Q 2 ''. This 
leads to the relation: 

rQ 2 

xg(x,Q 2 )= dk^^(x,k 2 T ) (12) 
Jo 

We can re-arrange eq.(|6|) in order to examine the azimuthal distribution of the final state 
partons: 

dyidy2d(p / aw J kx cos 9 — p T sin <p dt 

(13) 

The angle between px and p' T is and 

cos 9 = — sin 2 ± [ 1 — 7v sin 2 cos (14) 
k T \ k T j 

In the limit of k\ — > 0, the integrand becomes a delta function at the only allowed point in 
phase space, i.e. where = it. 

At this point we should emphasise that if one intends to use the standard DGLAP evolution 
equations to determine the Q 2 -evolution of xg then it is formally inconsistent to construct $ 
by taking the derivative and then using the &>r-factorisation formula above. In the DGLAP 
formalism, the QCD corrections to the cross section contain mass singularities, which must be 
removed by factorisation. This should be contrasted with the BFKL formalism where there 
are no mass singularities. However, we acknowledge that the absence of such singularities does 
not guarantee sensible results and if our cross section depends critically on inherently infra-red 
physics we may not be able to trust the results. 

These issues are dealt with nicely in the paper of Collins and Ellis [[[J. They introduce a 
factorisation scale to isolate the infra-red physics and, working in moment space, i.e. 

$(j kl; f/) = I dxx^^ix^a 2 ) (15) 
Jo 



they find that 

1 { k 2 \ 7c(j) 

$0; 4; /i 2 ) = 7cO') (^J xsf(j, /^) (16) 

where 7 c (j) is the BFKL anomalous dimension function which takes on the value of (exactly) 
1/2 when j — 1 + 12o; s ln2/7r. The BFKL corrections which multiply the gluon distribution 
are calculated by solving the BFKL equation with an infra-red cut-off of /i. To obtain the final 
cross section, $ must be convoluted with the 'impact factor' - in our case this leads to the 
expression: 

d2(jdir - f dv i~f, (~) xa ( x a 2, da{ ig ^ qq) 
d yi dy 2 ~ J dPT z Ji/e{z)xg{x,fi ) £ 

+ / dp% jf a dkl ^ j(u,£, 4, 9) zf l/e (z) *>(x, 4; Ml9 * d f m) (17) 

In this letter, we do not wish to focus on the specific solution for the function, $. Rather, 
we want to investigate the general features of this new factorisation prescription - and to 
ask whether it might be expected to produce significant deviations from the more standard 
formalism. 

To this end, let us use the following BFKL motivated form for the gluon distribution function 
(we ignore quarks throughout this paper - they are not included in the BFKL formalism since 
they do not lead to the dominant small x behaviour). 

=^>M- 2A-S«U ) (18) 

where 

A = ^41n2 

TT 

X" = ^28C(3) (19) 

7T 

The parameters, x and k^, are treated as essentially free - and we show results for different 
choices. Let us briefly reflect upon our prejudices regarding their likely values. The scale, /cq, is 
determined by the size of the hadron, and we therefore expect rather small values, ~ Aq CD . In 
the limit of x — >• 0, eq.(18) reduces to that which is predicted by the BFKL formalism. However, 
in order to ensure sensible behaviour at intermediate values of x, we modify the asymptotic 
logarithm, i.e. In 1/x — * \u{xq/x +1). This modification ensures a narrow distribution in k 2 at 
large x (consistent with the Q 2 -factorisation assumption that the partons are on-shell). One 
can think of xq as the parameter which delineates the onset of BFKL behaviour from the more 
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traditional large x behaviour, which suggests xq ~ 0.01. The normalisation, Af, is a constant 
and unimportant for the purposes of this paper. 

The following figures show the effect of /^-factorisation on the cross section for direct 
photoproduction of dijets at HERA, i.e. 820 GeV protons colliding with 30 GeV electrons. 
The cross sections are ep cross sections, and the photon flux was approximated using the same 
Weiszacker- Williams formula as in ref.fl]]. We took A = 0.5 and A" = 6.1, in eq.(19) but allowed 
a running coupling (with Aqcd = 200 MeV and Q 2 = in the hard scattering cross section. 

In figs.(l), we show the variation of the cross section with the jet rapidity (the jets are fixed 
to have equal rapidity, i.e. y\ = y 2 — y) for four different combinations of the parameters k 
and xq. In all cases, we find that the /^-factorisation prescription (eq.(17)) produces a cross 
section which is significantly smaller than that which is obtained using the Q 2 -factorisation 
prescription (eq.(3)) for those jets which are produced in the electron direction, i.e. small x. 
Not surprisingly, the effect increases as ko and Xq are increased - since we move the peak of the 
/^-distribution to larger energies as ko increases and broaden the width of the /^-distribution 
as xq increases. 

The transition from the k?- to Q 2 -factorisation prescriptions is illustrated in figs. (2). As the 
factorisation scale, n 2 , is increased the cross section tends toward the standard, Q 2 -factorisation, 
prescription result (shown as the horizontal dotted line on each plot). Exact agreement is 
never quite obtained since the upper limit on the /c-r-factorisation integral is calculated using 
the exact kinematics (rather than assuming it to be equal to Q 2 = which is the hard 

scale taken in the Q 2 -factorisation calculation). In fig. (2a), we show the /i 2 dependence at 
y = ±1 using the BFKL motivated form for $ and typical values of ko and xq. In fig. (2b), 
we show the equivalent plots using a $ which has been computed by taking the derivative 
of the GRV structure function parametrisations for the proton PD[ . As we remarked earlier, 
such an approach is strictly inconsistent with the BFKL formalism - one extracts a parton 
/^-distribution from the Q 2 -factorisation approach (in which the GRV structure functions are 
derived) - it is not clear that such an extraction is meaningful. Even so, our results again show 
significant deviation from the standard approach, especially in the backward direction. 

From the plots of figs.(l) and (2) it should be clear that the essential BFKL prediction that 
the incoming partons have non-negligible transverse momenta has important consequences at 
HERA. It destroys the notion of universal small x parton distributions in the standard (Q 2 - 
factorisation) formalism. Consequently, it is essential to measure the small-x parton distri- 
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bution functions in as many independent processes as is possible, e.g. xg(x,Q 2 ) from dijets 
JTJ], dF 2 /dlnQ 2 |TIJ, F L (T2| and J/\I/ production |T3| - we do not expect universal agreement 
between all processes. 

Another consequence of the non-zero parton kx is illustrated in fig. (3). The azimuthal 
distribution of the final state partons is plotted. Instead of back-to-back 'jets', we expect a de- 
correlating effect due to the imbalance of the incoming k?- This could show up experimentally, 
as a broadening of the azimuthal distribution for jets produced in the backward direction 
compared to those produced in more forward directions. We note that the apparent spikes 
at the peaks of the distributions are, of course, not physical and are simply due to numerical 
limitations. 

In summary, we have performed a quantitative study of the '/c^-factorisation' prescription 
when applied to a process that is currently being measured at HERA. Comparison with the 
traditional 'Q 2 -factorisation' approach shows that, in the region where small-x gluons within 
the proton are being probed, the discrepancy can be sizeable, i.e. ~ 20%. This deviation 
arises essentially due to the relaxation of the assumption that the partons which enter the hard 
scatter are on-shell and has the important consequence that we expect xg(x, Q 2 ) (for small x) 
to be process dependent. Another signal for such corrections would be a deviation from the 
ideal 'back-to-back' configuration in the azimuthal angle, as shown in fig. (3). 
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Figure Captions 

Fig. 1 A comparison of the Q 2 - and /cr-factorisation calculations of the dijet rapidity distribu- 
tions. We have chosen y± — y 2 — y and define negative rapidity to be in the proton 
direction. The normalisation of $ was adjusted to reproduce approximately the value of 
xg(x,Q 2 ) of GRV at x=0.01 and Q 2 = 10 GeV 2 . 

Fig. 2 To illustrate the transition from k?- to <5 2 -factorisation as the delineating scale, fi 2 , is 
increased (strictly speaking we use min(/z 2 , Q 2 ) as the delineating scale). In each case the 
prediction of eq.(3) is described by the dotted line, the first (second) term in eq.(17) by 
the long (short) dashed line and their sum by the dot-dashed line. 

Fig. 3 The azimuthal distribution of the final state partons, determined at y\ — yi — 1. and for 
two different values of k . 
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